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Justification of the choice of the subject: 


A suitable facility layout is with no doubt a key to efficient production 
that should be considered in the earlier stages of a project. A good 
layout might boost the efficiency and the effectiveness of the works. 
The objective is to arrange the temporary facilities such as the job 
office, the warehouse or the labor residence in a way that maximum 
design quality is achieved and the total cost of interactions among 
these facilities is minimized. An effective construction site layout 
planning is of the utmost importance for the success of any 
construction project. 
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l. Introduction: 


Suitable facility layout is believed to be the heart of efficient 
production that should be considered early in the planning phase. An 
appropriate construction site layout boosts the effectiveness and 
efficiency of the works. However, the arrangement of site facilities is 
hindered by many constraints such as limitations on the site area, 
adjacent buildings, access, the location and orientation of the 
building to be constructed. 


A Construction site layout problem can be modeled as a quadratic 
assignment problem (QAP) when the costs associated with flow 
between departments are assumed to be linear with respect to 
distance traveled and quantity of the flow. QAP problems are known 
as non-polynomial hard problems and because of the combinatorial 
complexity, these cannot be solved exhaustively for reasonably sized 
layout problems. As an example, for n facilities, the number of 
possible alternatives, that is the number of feasible configurations, is 
n! with larger growth than e”. Even for a small n. For 10 facilities, the 
number of possible alternatives is already well over 3,628,000. For 15 
facilities, we are already in the 12-digit numbers. In real problems, a 
project with n=15 can be considered as a small project. Due to the 
complexity of the site layout problems, numerous techniques have 
been proposed to find solutions to these problems; however, it is 
very difficult to obtain an optimal one suitable for hand calculations. 
Thus, optimization techniques seem to be suitable means to search 
for solutions of the site layout problems. 


In this research, two recently developed metaheuristic algorithms, 
known as colliding bodies optimization (CBO) and Enhanced colliding 
bodies optimization (ECBO), are applied to the solution of 
construction site layout problems. Colliding Bodies Optimization is 
developed by Kaveh and Mahdavi and Enhanced Colliding Bodies 
Optimization is presented by Kaveh and Ilchi Ghazaan. CBO and ECBO 


are employed for solution of the CSLP problem and results are 
compared with those of some previous algorithms. 


ll. Construction Site Layout Problem 


Construction site layout planning problems can be modeled as a 
QAP (Quadratic Assignment Problem) in which costs associated with 
the flow between facilities are linear with respect to the distance 
traveled and quantity of the flow. The objective of construction site 
layout planning is to assign a number of predetermined facilities (n) 
uniquely into a number of predetermined locations (m) where the 
number of locations should be equal or greater than number of 
facilities. If the number of predetermined locations (m) is greater 
than the number of predetermined facilities (n), then (m-n) facilities 
can be added to make both numbers equal. By assigning both the 
distance and frequency as 0, the facilities added will not affect the 
layout results. 
If each of the predetermined places is capable of accommodating any 
of the facilities, then the facility layout problem can be modeled as 
an equal-area facility layout problem. If some of the predetermined 
places are only able to accommodate some of the facilities, then the 
problem becomes an unequal-area facility layout problem, where 
predetermined places have differing areas. Generally, unequal-area 
layout problems are more difficult to solve than equal-area layout 
problems, primarily because unequal-area layout problems introduce 
additional constraints into the problem formulation. 


1. Objective function: 


The objective function of several models given in the Table below 
takes the general form: 
n 


n 
Minimize F = > 2 Wij * dij 
j=1 


i=1 


Where: 


e Fis the objective function and n is the number of facilities 
and locations. 


e Coefficient Wij represents either the actual transportation 
cost per unit distance between facilities i and j (taking into 
consideration the number of trips made) or a relative 
proximity weight that reflects the required closeness 
between facilities i and j. 


e dj is the distances between facilities i and j 


























No. | Pseudo model of the objective function 

1 To minimize the frequency of trips made by construction 
personnel 

2 To minimize the total transportation costs of resources 
between facilities 

3 To minimize the cost of facility construction and the 
interactive cost between facilities 

4 To minimize the total transportation costs of resources 
between facilities (presented through a system of proximity 
weights associated with an exponential scale) 

5 To minimize the total transportation costs of resources 
between facilities and the total relocation costs (presented 
through a system of proximity weights and relocation weights) 





Different kind of objective functions in the previous researches 
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. Layout representation: 


Each layout alternative can be represented by a n*n permutation 


matrix 


(A permutation Matrix is a square matrix that has exactly one 


entry of 1 in each row and each column and Os elsewhere), where n is 
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the number of facilities or locations, whose rows and columns 
represent facilities and locations, respectively. The permutation 
matrix allows a single entry of one in each row and each column, with 
all remaining entries being zero. The table below shows an example 
of a permutation matrix with 8 facilities and 8 locations. 


Each layout alternative can be represented by a n*n permutation 


L1 L2 L3 L4 L5 L6 L7 L8 L9 





Where Li is the it" location and Fj is the j'" facility. 


A specific solution to the site layout problem as shown by the table 
above is a very sparse matrix and would therefore consume 
unnecessary computing resources if it is used for large and practical 
problems. Therefore, because of the property of one-to-one 
correspondences between facilities and locations, a sequence of 
integers can be used as a more efficient alternative, like that: 


F1 F2 F3 F4 F5 F6 F7 F8 


Pelee eee 


Each position or entry in the sequence represents a facility; the 
integer number in the entry represents the location to place the 
corresponding facility. However, the sequence-based representation 
may lead to infeasible solutions where multiple entries in the 
sequence have the same integer number, specifically, the situation of 
overlay, when adopting the metaheuristic methods. Therefore, some 
modifications should be made to overcome this difficulty. 


lll. Metaheuristic Algorithms 


To solve construction site layout problems, we use two newly 
created metaheuristic algorithms, the colliding bodies optimization 
and the enhanced colliding bodies optimization. 


1. Colliding Bodies Optimization 


Colliding Bodies Optimization (CBO) is an efficient metaheuristic 
optimization algorithm that is based on one-dimensional collisions 
between bodies. All of the following explanations about this method, 
including definitions and formulas, are extracted from Kaveh and 
Mahdavi who presented this method. 

In this method, one object collides with the other object, and they 
move toward a minimum energy level. Collisions between these 
objects are governed by two laws of physics: momentum law and 
energy law. 

Each solution candidate Xi containing a number of variables (Xi={Xij}) 
is considered as a colliding body. The objects have assigned masses 
and are divided into two equal groups: stationary and moving 
objects, where the moving objects move to follow stationary objects 
and a collision occurs between pairs of objects. This takes place for 
two purposes: 


e To improve the locations of moving objects 


e To push stationary objects toward better locations 


After the collision, new locations of colliding bodies are updated 
based on the new velocities using collision laws. 


The Colliding Bodies algorithm can be briefly presented in the 
following 7 steps: 


Step 1: Initialization 


The algorithm starts with a random initial population of agents (CBs) 
in an m-dimensional search space by the following formula: 


x, = Xmin + rand°(Xmax = Xinin 
Where: 
e x; determines the initial value vector of the ith CB. 


© Xmax and Xmin are the minimum and the maximum allowable 
value vectors of variables. 


e rand is arandom number in the interval [0, 1]. 
e nis the number of colliding bodies. 
Step 2: Defining Mass 
Each colliding body (CB), Xi, has a specified mass defined as: 
n 1 
pO) 
fit(k) 
Where: 
e fit(k) represents the objective function value of the it” CB. 
e nis the number of colliding bodies. 
Step 3: Creating groups 


Objective function values are arranged in an ascending order. The 
sorted CBs are divided into two equal groups: 


e The lower half of the CBs are stationary CBs that have lower 
objective function values. These CBs are considered as good 
agents. 


e The CBs of the upper half are moving ones. These CBs move 
toward the lower ones and then the agents with upper value of 
each group collide together. 


Step 4: Criteria before the collision 


The initial velocities of stationary CBs are equal to: 


n 
vi = 0,i = 1,2, "5 
The velocities of moving CBs before collision are equal to: 
n n 
vi = ec = at Lot 2,..,n 


Where: 


e viand x; are the velocity and location vector of the ith CB in this 
group, respectively 


© Xi-n/2is the it CB pair location of x; in the previous group. 
Step 5: Criteria after the collision 


After the collision, the velocity of stationary CBs (V) are specified by: 


1 


(m,n + Em, n)v, n 
= 2 2 2 
vj) = ———— 


Mn + Mi 


where £ is the coefficient of restitution (COR) that decreases linearly 
from unity to zero. 


iter 


where iter and itermax are the current iteration number and the total 
number of iterations for optimization process, respectively. 
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Step 6: Updating the Colliding Bodies 


New locations of the CBs are evaluated using their velocities after the 
collision. The new locations of stationary CBs are: 


xit" = xi t+rand°v;’ 


Where: 


e xi" xi, and v are the new location, previous location, and the 
velocity after the collision of the i CB, respectively. 


e rand is a random vector uniformly distributed in the range of 
[1,1]. 


e The sign ° denotes an element-by-element multiplication. 
Step 7: Termination criterion check 


The process of CBO algorithm is repeated from Step 2 until a 
termination criterion, such as maximum iteration number, is 
satisfied. 


2. Enhanced Colliding Bodies Optimization: 


Enhanced Colliding Bodies Optimization (ECBO) is a new version of 
the CBO which improves the CBO to get faster and to obtain more 
reliable solutions. Unlike CBO, the main feature of the ECBO is that it 
uses a memory to save some best solutions that cause an increase in 
the convergence speed of ECBO with respect to standard CBO. In 
order to improve the exploration capabilities of the CBO and to 
prevent premature convergence, ECBO utilizes a mechanism to 
escape from local optimal. 

All of the following explanations about this method, including 
definitions and formulas, are extracted from the works of A. Kaveh 
and |. Ghazaan. 

In order to introduce the ECBO, the following steps are developed: 


Step 1: Initialization 
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The algorithm starts with a random initial population of agents (CBs) 
in an m-dimensional search space by the following formula: 


x; = Xmin + rand? (Xmax — Xmin) 
Step 2: Defining mass 


The value of mass for each CB is evaluated according to the same 
equation as before. 


Step 3: Saving 


In this step, colliding memory (CM) is utilized to save a number of 
historically best CB vectors and their related mass and objective 
function values with the aim of improving the algorithm’s 
performance. At each iteration, solution vectors that are saved in the 
CM are added to the population and the same number of the current 
worst CBs are deleted. Finally, CBs are sorted according to their 
masses in a decreasing order. 


Step 4: Creating groups 


The CBs are divided into two equal groups according to their 
objective function values: stationary and moving group. 


Step 5: Criteria before the collision 


The velocities of the stationary and moving bodies before collision 
are evaluated by the same equations used in the CBO. 


Step 6: Criteria after the collision 


The velocities of the stationary and moving bodies after collision are 
evaluated by the same equations used in the CBO. 


Step 7: Updating CBs 


The new location of each CB is evaluated by Equation used in step 6 
of the CBO. 


Step 8: Escape from local optimal 
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In order to escape from local optimal, a parameter like Pro within (0, 
1) is introduced, which specifies whether a component of each CB 
must be changed or not. For each colliding body Pro is compared 
with rnj(i=1,2,...,n) which is a random number uniformly distributed 
within (0, 1). If rni < pro, one dimension of the it CB is selected 
randomly and its value is regenerated as follows: 


Xij = Xjmin + random * (Xj max — Xjmin) 
Step 9: Termination criterion check 


After the predefined maximum iteration number, the optimization 
process is terminated. If this criterion is not satisfied, go to Step 2 for 
a new round of iteration. 


IV. Discussion of the results 


In CBO and ECBO algorithms, each solution candidate Xi 
containing a number of variables is considered as a colliding body. 
Each CB is considered as a sequence of variables that represents a 
layout solution and different sequences mean different layout 
solutions. Each variable in the sequence represents a facility, and the 
value of variable indicates the location that is assigned to the 
corresponding facility. Since every location is capable of receiving 
only one facility, the CBs should not have duplicated values; violation 
from this point generates infeasible solutions. However, all the 
variables of a CB in CBO and ECBO are independent of each other; 
thus, updating the velocity and position of a CB is performed 
independently. Therefore, more than one variable in an updated CB 
may have the same value. Thus, some modifications in updating 
mechanism should be performed to overcome this infeasibility. The 
updating mechanism of the CBs is explained in the following: 
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Sort variables’ velocity of CB 







N€ Number of variables 
Cé 1 





Has the updated 
value of CV been selected 
before? 


Revert updated value of 


CV to previous value 


V. Case Study 


To show the applicability and performance of the CBO and ECBO 
algorithms for construction site layout optimization a case study is 
shown below. Their results are compared of those of the PSO 
(Particle Swarm Optimization — a classic optimization method) The 
algorithms are coded in MATLAB R2011a, and the experiments are 
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performed on a personal computer with Intel®CoreTM i7 processor 
(1.73 GHz) and 4 GB RAM under the windows 10 Home 64-bit 
operating system. Parameter values used in these case studies are 


shown in the table below: 





PSO CBO 





Population size: 50 Population size: 50 
Inertia Weight:0,4-0,9 
C1=C2: 2 











ECBO 
Population size: 50 
CM Size: 5 
Pro: 0,3 








Results are as follows: 


This case study is a medium-sized project and is taken from Li and 
Love. The purpose of this problem is to find the most appropriate 
arrangement for placing 11 facilities into 11 predetermined locations 
on the site. The 11 facilities and their corresponding index numbers. 


In this case study are shown in the table below: 















































Index Site facilities Note 

number 

1 Site office Not fixed 
2 False work workshop Not fixed 
3 Labor residence Not fixed 
4 Storeroom 1 Not fixed 
5 Storeroom 2 Not fixed 
6 Carpentry workshop Not fixed 
7 Reinforcement steel workshop Not fixed 
8 Side gate Fixed to 1 
9 Electrical, water, and other utilities control room Not fixed 
10 Concrete batch workshop Not fixed 
11 Main gate Fixed to 10 





For the construction site layout selection, two assumptions are made: 
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e Each of the predetermined locations is capable of 
accommodating any of the facilities. 

e The main gate and side gate are treated as special facilities, 
which have been fixed on the predetermined locations. 


Objective function: 


The objective of this case is to minimizing the total traveling distance 


of site personnel between facilities: 
n n 
` ` Xin Xfi dei 


1 l=1 


M= 


n 
minimize TD = ` 
i=1 


= 


<=. 
II 
m. 


3 


With: i Xij = i Aati = 1 


~ 


Where: 


e nis the number of facilities, in this case n=11. 

e xi=1 if the facilty | is assigned to the location k and equals O 
othrrwise. 

e xis defined similarly. 

e Coefficient fj is the frequency of trips made by construction 
personnel between facilities i and j per day. 

e Coefficient dx is the distances between the locations k and |. 


Distances between site locations: 


tocaion| af taf at | sf tat tot 
1l 





Trip frequencies between the facilities: 





Results: 


This example is solved by carrying out 50 independent optimization 
runs through 200 iterations to obtain statistically significant results 
by PSO, CBO, and ECBO. Results returned are: 


Algorithm Standard deviation 
12,546] 12,56| 12,756 47,39 





As it can be seen from this table, the average, worst, and standard 
deviation for ECBO are: 12,555, 12,746, and 32.11, respectively, 
which are better than those of CBO and PSO. This indicates that ECBO 
not only finds a better best solution, but also it is more stable. The 
convergence curves for the ECBO, CBO, and PSO in terms of the 
number of iterations are shown in the figure below: 
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VI. Conclusion 


In the study above the application of two recently developed 
metaheuristic algorithms, CBO and ECBO, is introduced to solve 
construction site layout problem. The governing laws of physics 
initiate the base of the CBO and ECBO algorithms, where these laws 
determine the movement process of the objects. CBO utilizes simple 
formulation to find minimum of objective functions and does not 
depend on any internal parameter. In order to improve the 
exploration capabilities of the CBO and to prevent a premature 
convergence, ECBO uses a mechanism to escape from local optimal. 
The latter also uses a Colliding Memory to save a number of the so 
far best solutions to reduce the computational cost. To validate the 
models, a case study is considered. The result verify that the 
proposed approach performs very well both in finding better results 
and using lower number of evaluations to find the optimum. 
Comparison of the results with some other well-known 
metaheuristics shows the suitability and efficiency of the utilized 
algorithms in CSLP. The proposed algorithms are highly competitive 
with other metaheuristic algorithms in quality of solutions and 
convergence speed. 
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VII. Implementation in Morocco: 


After all that a question might come to mind, is this optimization 
method applicable in Morocco? 


A brief and simple answer is yes, it is applicable and its application is 
quite easy and simple. 


In fact, using the 2 metaheuristic Algorithms presented in this 
research; the CBO and the ECBO isn’t a complicated task neither an 
expensive one, the proposed algorithms work in almost the same 
way as those already used in Morocco so engineers and technicians 
must have enough knowledge to make it work, and like said before it 
will boost efficiency and effectiveness and minimize the cost, and 
with these being the main objectives of the work of an engineer | 
think it is worth it to apply this in Moroccan construction sites. 
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